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Abstract 

It is shown that nonzero orbital momentum in the vertex of secondary inter- 
action in the triangle graph leads to a more clear picture corresponding to the 
moving complex singularity compared with the case of constant vertex. The peak 
in the distribution of the final particles invariant mass becomes narrower and its 
shift with momentum transfer more distinct. It opens the possibility to observe 
this picture experimentally. 
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The triangle mechanism (see Fig.l graph) is one of the simplest mechanisms of direct 
nuclear interactions. It is very frequently used in the description of nuclear reactions. A 
characteristic graph of this type corresponds to a primary interaction (the right upper vertex) 
followed by secondary elastic or inelastic interaction (the lower vertex). The experimental 
determination of its contribution as well as separation of the kinematical regions of its 
domination would permit us to extract valuable and reliable information on nuclear structure 
and the off-shell amplitude of the "elementary" process in the lower vertex. 

Until now very few cases are known when the dominant contribution of the triangle 
mechanism is seen. Let us point to the peak in pion spectrum from the capture of stopped 
kaons on the deuteron K~d — > pA7r _ [1] due to the triangle graph with NE intermediate state 
and the conversion NE — > pA [2] . Besides it is possible to mention anomalies in photonuclear 
interactions [3] and enhancements in two nucleon mass spectra near 2mj^ due to FSI [4]. 
On the other side, it is well known that the characteristic feature of the triangle graph in 
nuclear physics consists in so called moving complex triangle singularity. Its position in 
the invariant mass M of the lower group of final particles depends on the value q of three- 
momentum transfer from initial to final fast particles in the right vertex [5] . The distribution 
with respect to M depends on the value of q. It is a model-independent criterion of the 
dominant role of the graph [6]. Up to now this test was not used for the identification of 
the triangle mechanism. The estimations for the process pd — > pdr] [7] and for several other 
processes with an intermediate A-isobar [8] predict a clear picture of the moving singularity 
but there will be rather broad maxima in the M distributions with widths of the order of 
100 -T- 200 MeV. 

Of course, the most distinct picture could be expected in the region of the excitation 
energy E ex = M — (mi + m 2 ) of the order of nuclear binding energies and momentum 
transfer q of the order of nuclear Fermi-momenta. Here the complex triangle singularity 
comes close to the physical region and narrow peaks and background suppression could be 
anticipated. Moreover, the cross sections in this region are much larger. It is hindered by 



the fact that Fig.l triangle graph has also the threshold root singularity at E ex = and 
usually it determines a picture of the spectrum at small £ex- As we shall see below, the 
picture of the moving singularity get too smooth and hardly yields to observation. 

However, the situation sharply varies if not the s-wave but p- or higher waves dominate 
in the amplitude of secondary interaction. The shape of the cusp singularity changes and 
it ceases to dominate at small E ex - As will be shown below, the situation with observation 
of the moving triangle singularity is essentially improved and we can see very clear and 
unambiguous picture. Later on it will be illustrated on graphic examples. Yet at first we 
will develop a formalism for Fig.l graph calculation with arbitrary nuclear form factor in 
1 + 3 vertex and account of nonzero orbital momenta in the lower vertex. 

Except the quantities M, E ex and q = p x — p y (all momenta in lab. system), introduced 
earlier, we will also use the nuclear binding energy e = mi + m 3 — and the characteristic 
momentum k in the nuclear vertex A — > 1 + 3: k 2 = 2mi 3 e, where TO13 is the reduced mass 
of particles 1 and 3. We will usually consider nuclear reactions and take the particles A, 1 
and 3 to be nonrelativistic. It is convenient to use the dimensionless variables £ and A [9,7], 
associated with M and q: 

c-c(m \ m A \ m iq 2 M 2 + m 2 -/j, 2 \ 

A = A(M, q) = (2) 

These expressions are valid in general case of relativistic particle 2. In the nonrelativistic 
case £ is expresses only through M: 

^(M) = ^^ M - mi ~ m2 . (3) 
m 3 M e w 

In terms of these variables the triangle graph has a square root branching point at £ = 
and logarithmic singularities at £ = A — 1 ± iy/X. 

Let us go to the estimation of Fig.l graph for the case of arbitrary upper left (nuclear) 
and lower (secondary interaction) vertices. We will take the upper right (primary inter- 
action) vertex to be constant, equivalent to factoring out its value in the triangle integral 



at some point near the maximum of the integrand. Usually this means that we factor out 
the amplitude of the process x + 3^y + 2at the proper momentum transfer and at the 
energy corresponding to the interaction of particle x with particle 3 at rest. The amplitude 
of the graph has the form of a four-fold integral over the 4-momentum p of particle 3 [10]. 
The integration over the fourth component reduces to determining the residue in the pole 



of the propagator of particle 1. (Then neglected terms are of the order y e/m or e/ fi where 
m and /i are nucleon and pion masses.) So in reality we have three-fold integral over p, 
containing propagators of particles 2 and 3 and two nontrivial vertices. The nuclear vertex 
A — > 1 + 3 together with the propagator of particle 3 corresponds to a nuclear wave function 
^iip)Yim(p)- The lower vertex can be represented as 

Mk)Yt M {k) (4) 

where k is the relative momentum of particles 1 and 2 

k = p + Q, Q= 7 q. (5) 

mi + m 2 

Here I and L are the orbital momenta of the nuclear vertex and the secondary interaction 
vertex. As for the latter, we deal here only with nontrivial part (4) of its amplitude which 
contains the dependence on k. The most interesting phenomena appear in the M region not 
far from the threshold of the channel (1+2), i.e. for mi + m 2 not far from M, where the 
dependence on k is most important. Actually all vertices also contain spin parts (Clebsch - 
Jordan coefficients due, for example, to spin-orbital interaction), but these do not influence 
the integration procedure. 

The propagator of particle 2 can be transformed to the form 

2m 2 E 2 -pl + ir) = _ mi + m2 {( p + Q) 2 - 2m 12 E ex - i V } . (6) 

mi 

Note that Q 2 = A/t 2 , 2m\ 2 E ex = £k 2 . The propagator (6) contains the same vector combi- 
nation as in eq.(5) that is the relative momentum of particles 1 and 2. This is important for 
our formal transformations. This statement is valid also in the case of a relativistic particle 2 



(see ref. [9]). The expression for the amplitude of the triangle graph of Fig.l is proportional 
to 

T _ f , %ip)YUP)M\P + Q\)Y£ M (p + Q) m 
A ~J P (P + Q) 2 - 2m 12 E ex - if] ■ [n 

Here E ex can take positive as well as negative values. 

Let us transform eq.(7) to coordinate space and use nuclear wave function in coordinate 
space ipi(r)Y lm (r): 

%(p)YUp) = J Mr')Y lm (i')e ipr 'dr'. (8) 

If the lower vertex were constant we would use one of the following relations to transform 
the denominator of eq.(7): 



dr ' TT~lS5 = ir\ dr - ( 9 ) 



k 2 — a 2 — irj Att J r k 2 + (3 2 An 

Then simple calculation would lead to the result (for example if E ex < 0) 

roc I 

J A = (2tt)¥>U-Q) / drr 3l (VXKr)exp(-^-^r)^(r) . (10) 

Jo v 

For the case E ex > the quantity — must be substituted for 

For an arbitrary lower vertex we will apply a transformation to coordinate space which 
includes a new function fL^r) which will be described later: 

<p L (k)Yl M (k) 



k 2 + (3 2 

Then straightforward but tedious calculations lead to the result 



jYl M (r)f L (r)e-^dr. (11) 



I A = (270=7^(2/ + l)(-l)^^QoL Q^ t (-Q) (12) 

T,t 

roo 

■ / drr 2 j T (Qr)i)i(r)f L (r) 
Jo 

where jt is a spherical Bessel function. The nontrivial part of I a, which will be denoted as 
M 

oo 

M= jdr-r 2 j T (Qr)^(r)f L (r), (13) 



contains one-fold integral with radial part of nuclear wave function and function fi(r) asso- 
ciated with the lower vertex of the triangle graph of Fig.l. Here T, / and L must obey the 
triangle rule. The relative weights of terms with different T in the final results for the cross 
section, polarizations etc. will depend not only on factors in eq.(12) but also on the spin 
structure of all vertices in the triangle graph. 

Let us turn to eq.(ll). Inverse transformation gives 

Mr) = £]zjggle«. (U) 



Various choices of this function are described in ref. [9]. Here we consider only one kind 
chosen such as to give a correct behavior of the amplitude for the case of arbitrary orbital 
momenta in the near-threshold region: 

<p L (k)=Ck L . (15) 

Then for example 

AM = ~- = -^Af'WH (16) 

In the general case, a minimum account of the centrifugal barrier requires 

fL(r) = (-l)WC^hP(iPr), (17) 

where is the spherical Hankel function of the first class. Expressions (11), (14), (16), 
(17) correspond to £ < when (3 = >/— For the case £ > the quantity (3 in these 
expressions must be changed into —27 where 7 = y/^K. 

The formula (13) enables us to carry out numerical estimations in a simple way. Here we 
will show the results of concrete calculations for several cases with s- and p- orbital momenta. 
For s-wave nuclear wave functions we used the deuteron function in Paris potential as well as 
the wave function in Gauss parametrization ip(r) ~ exp(— p\r 2 j2\ For the p-wave we used 
a p-wave Gauss function ip(r) ~ rexp{— p 2 ) r 2 /2) as well as a "quasi-Hulten" p-wave function 
of the form 

Mr) ~ ( — + -L.) ( e ~ Kr - 3e-^ r + 3e^ K+2 ^ r - e ^ K+ ^ r ) (18) 
\Kr K z r z J v ' 



which has correct asymptotic behavior at r — > and r — > oo. 

Fig. 2 shows the results for |.M| 2 for the case of the deuteron with Paris wave function: 
a) for s-wave lower vertex and b) for p-wave lower vertex. In the latter case there is now 
a distinctive and pronounced systematic pattern in which the maximum of the distribution 
moves to larger values of £ with increasing A. For s-wave Gauss function the pattern are 
close to the deuteron case. Fig. 3 illustrates the corresponding effect for the case of p-wave 
Gauss nuclear function with p = 0.75k (it corresponds to 12 C). Again the picture is more 
distinct for the case of p-wave lower vertex. The p-wave "quasi-Hulten" function (18) gives 
similar results. These calculations were done with the amplitude of lower vertex of the Fig.l 
graph in the form (22). Introduction of cut-off factors (see ref. [9]) does not lead to the 
qualitative changes. We see that the cases with the p-wave lower vertex open the possibility 
to observe the picture of the moving triangle singularity experimentally. Up to now this 
picture has not been exploited for the identification of the triangle mechanism. 

In conclusion it can be said that we have first derived the simple and transparent formulas 
(12)-(14) for the triangle graph with arbitrary orbital momenta in the nuclear vertex and 
in the vertex of secondary interaction. The amplitude is represented in the form of one- 
fold integral in coordinate space. It is well suitable for numerical calculations due to the 
exponential decrease of the nuclear wave function. 

We then have pointed out that nonzero orbital momenta in the vertices have an important 
effect on the shape of experimental distributions. This idea is supported by the graphical 
results. The corresponding singularity manifests itself much more distinctly and sharply than 
for the case with L — 0. The experimental demonstration of this effect requires favourable 
situations. This is not easy, since the most interesting region is located close to the threshold 
for particles 1 + 2, where the p-wave amplitude is normally suppressed by the barrier. One 
interesting example is iViV-interaction in the final state for which the p-wave is greatly 
enhanced at low energies and can be separated by the choice of particular channels [11]. A 
second example is the 7riV-interaction which is dominated by p-wave part from rather small 



energies. 
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Figure 1: Typical triangle graph. 




Figure 2: \M.\ 2 as a function of £ for A = 1 (solid line), 3 (dotted line), 6 (dashed line), 9 
(dash-dotted line) in the case of the deuteron with Paris wave function: a) for s-wave lower 
vertex, b) for p-wave lower vertex. 




Figure 3: The same as in Fig.2 for p-wave nuclear function in the Gauss parametrization: 
a) for s-wave lower vertex, b) for p-wave lower vertex. 



